function f = elliptic_uy(x,y)
global alpha theta
f=0.2e1*(x.^2+y.^2+eps).^(-1+1/ alpha).*(y.*sin(0.2e1*(atan2(y+eps,x+eps)+theta)/alpha)+cos(0.2e1*(atan2(y+eps,x+eps)+theta)/alpha).*x)/alpha;